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Abstract 

In heavy-ion cohisions at very high energy the density of produced jets can be so high that 
the possibility of hadrons produced by recombination of shower partons in overlapping minijets 
may become important. We study such multi-minijet contribution to the hadron spectra and to 
dihadron correlation in Pb-Pb collisions at 2.76 TeV at LHC. We adjust the parameter controlling 
the momentum degradation of semihard partons by fitting the charged-particle distribution up to 
Pt ~ 16 GeV/c. The relative magnitudes of different identified hadrons and of various partonic 
components are fixed by the recombination formalism. We find that the coalescence of shower 
patons from adjacent miniijets can be as much as from single jets for meson production, and even 
more so for proton, but never dominant over other components. In 3-shower-parton recombination 
the ratio of 2-jet to 1-jet contributions increases with collision energy; its maximum can exceed 2 
at 5.5 TeV. Two-hadron correlation exhibits a broad peak on transverse rapidities, confirming that 
minijets play a central role at low px- 
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I. INTRODUCTION 



In heavy-ion collisions above 2 TeV the density of minijets produced by semihard scat- 
terings of partons can be so high that conventional treatment of such collisions may be 
inadequate. What is conventional at lower energy (0.2 TeV) is hydro dynamical description 
for transverse momentum pt < 2 GeV/c [l| and jet fragmentation at > 6 GeV/c [2]. In 
the intermediate region neither approach is reliable. The unconventional treatment based on 
quark recombination has shown some success in filling the gap js-S]. Conceptually, one can 
imagine that at extremely high energies, higher than what is available at the CERN Large 
Hadron Collider (LHC), heavy-ion collisions can be so explosive with a preponderance of 
minijets that hydrodynamics would be invalid even at low px on the one hand, and absence 
of jet-jet interaction would be unreasonable at moderate px on the other. The question 
is what the reality is at LHC where charged particle distribution is already known up to 
Pt ~ 18 GeV/c [6]. An answer cannot be given without some model-dependent analysis. 
Our aim in this paper is to find the degree of importance of multijet recombination in the 
formation of hadrons. 

Recently, the spectra of identified hadrons in Pb-Pb collisions at LHC for p^ < 5 GeV/c 
have been studied [?, 8|; it is found that the recombination of thermal and shower partons is 
important down to p^ ^ 1 GeV/c. It means that there is an abundance of minijet produced 
at ^/SNN■ = 2.76 GeV. The basic parameters (T, inverse slope of thermal partons, and k, a 
measure of the momentum degradation of semihard partons) are determined by fitting the 
data on vr, K,p and A distributions simultaneously. With those parameters at hand, we can 
then extend the study to higher px with sufficient constraint to confront the data on charged 
hadrons. 

Low-energy jets with transverse energy Et ^ 10 GeV are not the sort of objects that can 
be identified by jet-finding algorithms, such as HIJA jol, for heavy-ion collisions because 
they merge into the backgrounds from high-multiplicity underlying events in search for jets 
with Et > 50 GeV 10|]. Separation of neighboring minijets is not only unfeasible, but 
probably not meaningful. A discussion of jets usually starts with a definition of what the 
minimum energy in a jet cone (E^^^) is after background subtraction, and what the radius 
of the jet cone (R) is. It is based on the concept that particles in a jet are correlated 
with the initiating hard parton, and can be isolated from the background by an effective 
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algorithm. The task for accomphshing that becomes more and more difficult as the jet 
energy is reduced because of the fluctuation of the background. From the point of view of 
Pt distribution that is averaged over all pseudorapidity rj and azimuthal angle 0, at least 
in a midrapidity interval, say \ri\ < 1, and over all events, the background is the dominant 
part that is exponential. Small deviation from that part at a slightly larger px reveals the 
effect of semihard scatterings that cannot be easily identified on the event-by-event basis as 
well-defined jets in terms of Et and R. Nevertheless, we need to focus on them in order 
to calculate their effect on the pt distribution. For that purpose we use the term minijet 
without being precise about and R, since we shall not examine the event structure in 
the ■)]-(/) plot. We do consider semihard partons with transverse momenta A; > 3 GeV/c and 
the shower partons that they generate after emerging from the medium surface. Each such 
semihard parton and the cluster of associated shower partons will be referred to as a minijet. 

In heavy-ion collisions there are various theoretical issues related to minijets that have not 
yet evolved to a mature subject with general acceptance. The medium effects on semihard 
partons are important but hard to make precise, and the hadronization process is still con- 
troversial. The shower partons not only depend on the degree of momentum degradation in 
the medium, but also have various channels of hadronization, such as through recombination 
with thermal partons on the one hand and with other shower partons on the other. At LHC 
the high density of jets creates the possibility of shower partons from different jets overlap- 
ping in common spatial proximity so that their coalescence cannot be ignored. The study of 
multi-minijet contribution to the hadronic pt distribution and to two-particle correlation is 
the main concern of this paper in addition to the usual components involving single jets. We 
shall examine all possible components so as to exhibit the relative importance of each up to 
Pt ~ 16 GeV/c. From ^/s = 2.76 TeV with known charged-particle spectra, we extrapolate 
to 5.5 TeV to show how the multi-minijet contributions depend on collision energy. 

Although we can calculate the minijet contribution to the hadronic spectra, it is difficult 
to conclude from the observed pt spectra that minijets are necessarily existent. To make 
that conclusion cogent, two-particle correlation exhibited in terms of transverse rapidities has 
been used experimentally to show the existence of a broad peak [11!]. We shall analyze our 
calculated results in those variables and show general agreement with the data. Thus there is 
common ground in recognizing the important role that minijets play in heavy-ion collisions. 
However, our approach differs from the experimental approach in that we provide the par- 
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tonic basis of the observed phenomenon. At low pt where correlation is strongest, thermal 
partons play an important role despite the fact that they are by themselves uncorrelated. 
That is because two shower partons from one jet are correlated by momentum constraints, 
so their separate recombination with different thermal partons results in hadronic correla- 
tion. Thus Pt correlation among particles that are not too far apart in pseudorapidity rj 
and azimuthal angle is a fertile ground to find the footprints of recombination. 

Our study here is limited to midrapidity in central collisions and all our formulas are 
averaged over the azimuthal angles. As will become self-evident, the theoretical description 
that includes all components of hadronization (e.g., 7 components for proton production) is 
sufficiently complicated that it becomes essential to build first a clear and solid foundation for 
the physical processes without involving the azimuthal complexities. We shall use schematic 
diagrams to help the visualization of the various processes. 

In the next two sections we describe the multijet recombination processes for pions and 
protons. Two-particle correlations are discussed in Sec. IV. The parameters we use to do 
detailed calculations are given in Sec. V with results on single-particle spectra shown in Sec. 
VI. Extension to 5.5 TeV is made in Sec. VII. The results on two-particle correlation are 
presented in Sec. VIII. Conclusion is then given in the final section. 



II. TWO-JET RECOMBINATION FOR PION PRODUCTION 

All basic equations for the inclusive px distributions of hadrons (tt, K, p, A) produced in 
Pb-Pb collisions at LHC have already been given explicitly in Ref. |8|, whose equations will 
be referred to hereafter with^ the prefix I. The definition of all the quantities undefined in 
this paper can be found in [8]. We shall not repeat them here. Instead, we show schematic 
diagrams that can clearly depict the processes involved. They include both the types of 
processes already considered in js] and new ones, such as two-jet recombination, for which 
relevant equations will be given. 

In Fig. 1 we show the diagrams in the transverse plane for the recombination of (a) 
thermal T and shower S partons, (b) SS in one jet, and (c) SS in two jets, which will be 
abbreviated by TS, (SS)"'^-' and (SS)^-', respectively. In the notation of Eq. I-(35), k is the 
momentum of the semihard parton at creation, and q is the momentum at the medium 
surface. The thick red vectors have the dual role of representing the jet momentum in the 
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FIG. 1: (Color online) Schematic diagrams for parton recombination of (a) TS, (b) SS in one jet, 
and (c) SS in two jets. Thick (red) lines represent partons in medium, thin (red) lines partons out 
of medium, thinnest (red) lines shower partons, and dashed (blue) lines thermal partons. 



medium and the degradation effect described by G{k, q, ^) in the same equation 8|, ll2|. The 
thinner red lines outside the medium are the semihard partons qj, which can emit shower 
partons represented by the thinnest red lines denoted by pj. The blue dashed arrows are 
thermal partons. Recombination is presented by a large black blob with the outgoing open 
arrow depicting the produced pion. The lengths and angles of the vectors are not drawn to 
scale due to the limitation in presenting the figures clearly, and should not be taken literally. 

Figure 1(a) is described by I-(14), and Fig. 1(b) by I-(7), (18). For Fig. 1(c) to occur, 
the two vectors ki and ^2 should be nearby and approximately parallel so that there is an 
appreciable overlap of the jet cones. Considering only partons that are near 77 = 0, and 
exhibiting the (p dependence, we have for the two-jet contribution to the pion spectrum 



2 r 



Pj 



Rr(pi,<^i,P2,02,PT,0), (1) 



dpTt 

where 5(p, ^) is the integrated shower parton distribution defined in I-(6) in terms of the 
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semihard parton distribution Fi{q,^) at the surface and the shower distribution Si{p/q) in 
the i-jet, integrated over q and summed over i. The variable ^ is a measure of the dynamical 
path length in the medium created by the heavy-ion collision at impact parameter b 12|. 
P(^, (p, b) is the probability for ^ to occur for a path at angle (p initiated at {xq, yo), weighted 
by the nuclear overlap function, and integrated over all (xo,?/o)- The quantity in the square 
brackets in Eq. ([1]) is the probability of having shower parton pj in Fig. 1(c). is the 
recombination function (RF) characterized by F that summarizes all other dependencies 
besides pj and such as the spatial separation of the shower partons. It is important to 
recognize that the product (j = 1, 2) in Eq. ([T]) implies the creation and degradation of two 
independent semihard partons, but whose shower partons must be in close proximity outside 
the medium if they are to coalesce to form the pion. Thus only at high energy (as at LHC) 
is the jet density high enough for the 2-j recombination to occur. 

For the RF, Rp(pi, 0i,p2, 02,Pt, 0), the coalescence process clearly cannot take place if 
01 and 02 are not nearly equal, since non-parallel partons have large relative momentum 
transverse to pi +p2- Large relative longitudinal momentum parallel to pi +p2 is permitted 
in the parton model, since the momentum fraction of a parton in a hadron can vary from 
to 1. Relative momentum transverse to that is limited by the confinement restriction that 
it should not exceed the binding energy of the constituents. One may consider a Gaussian 
distribution in |0i — 02| with an appropriate width. However, since 0i and 02 are integrated 
over in Eq. ([T]), it is simpler to adopt a factorizable form that requires the partons to be 
parallel but with a suitable normalization factor F that we can estimate, i.e., 

R?(P1,01,P2,02,PT,0) = r5(0i - 02)5 (^^i±ia _ 0^ i?-(pi,p2,PT), (2) 

where F is the probability that two parallel partons can recombine. Since the partons are 
emitted from the medium at early times, we may consider the emitting system as being 
a thin almond-shaped overlap region viewed from its side in the same transverse plane at 
midrapidity as where the pion is detected. For centrality c < 0.05 the almond is almost 
circular. The partons at 0j are parallel, but can be emitted at any distance from the center 
of the circle. Looking at the emitting source edgewise, it is essentially a one-dimensional 
system of width approximately 10 fm, which is slightly less than 2Ra since high-density 
partons are not likely to be emitted tangentially from the edges. The two parallel partons 
should be separated by a distance not greater than the diameter of a hadron (~ 2 fm), given 
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PrdpT Pt J Pi P2 



— Fi{qa, K)Si{pa/qo 



RiPl,P2,PT), (3) 



that the jets have some width. Thus our estimate for F is the ratio ~ 2/10. We do not see 
that any more elaborate analysis of the coalescence process can provide a more transparent 
description of Rf. 

Applying Eq. ([2]) to ([T]), we obtain 

j-j dqa 

a=l 

where the integral over C,j in ([1]) has been replaced by the average distribution Fi{q,K), 
discussed in sj]. Further consideration of the parameter k will be explained below, but for 
now we give the original definition 

l/K = q/k (4) 

which is the fraction of the semihard parton's momentum k that is not lost to the medium as 
it emerges with momentum q. All parts of Eq. are known from [8] , so it is straightforward 
to calculate the 2-j contribution to pion production. 

It should be noted that theoretically it is possible for a very hard parton to split into two 
softer partons which can separately form two jets that are close-by and may even overlap. 
Such a QCD process would have only one Fi function at a much higher parton momentum. 
However, since the parton distribution fi{k) at creation decreases rapidly with increasing 
k, the production of such hard partons is suppressed compared to the semihard partons 
that we consider, as our results will show the dominance of minijets. For that reason such 
processes are not included in our calculation of dN^^ /pxdpT- 



III. MULTIJET RECOMBINATION FOR PROTON PRODUCTION 

There are many components of thermal and shower recombination that contribute to the 
production of protons. Symbolically, they can be expressed as in I-(19) 

F,,, = rrr + tts + t{ss) + {sss) + t{ss) + [{ss) ''s] + {sss) (5) 

Except for the first term that does not involve any S, the other six terms are depicted 
by the six figures in Fig. 2, respectively. The equations for TTS, T(SS)^'', and (SSS)"*^-^, 
corresponding to Fig. 2(a), (b) and (c), are given in I-(A9), (AlO) and (All), respectively. 
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T(SS)i^' 





T(SS)2J' 
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[(SS)i^S]2j' 



(SSS)3J' 



FIG. 2: (Color online) Diagrams for proton production by recombination of partons with same 
line- types as in Fig. 1. 



For the next three processes they are the following 



p 2Q+/3+3 

m^Pj. JO 



PT~Pl 

dpi I dp2{piP2T{pT-pi-P2Y 


\Y.P,e-^^'^S\pu.K)S\p,,K). 



X 



(6) 



{jM) 



dK 



[{SS)'-^S] 



1JC123 



Prdpi 



p 2a+/3+3 

■mij,prp ^ Jo 



dpi 



PT-Pl 



dp2{piP2T{PT -P1-P2Y 

^-Y,S'^\p,,Pk.K)S\p,,K). 

{jM) 



p 2a+/3+3 

m^Pj, Jo 



dpi 



PT-Pl 



dp2{piP2T{,PT - Pi- P2YW S^i^p 



(7) 



(8) 



PrdpT 

where ^(j^^) denotes cyclic permutation of {j,k,i) over (1,2,3) with ps = px — Pi — P2- 



The exponents a = 1.75 and (3 = 1.05 are from the proton RF 



The quantities 5^ 
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and iS^'' are defined in I-(A4) and (A12) for one quark and two quarks in a gluon jet. Tlie 
sum over all semihard partons is approximated by a times gluon contribution only with 
cr = 1.2 signifying that all other quark jets are regarded as contributing ~ 20% more to 
the gluon jet. This approximation is based on concrete calculations of certain components 
in which we compare the u quark contribution to that of the gluon. We have found that 
using a = 1.2 as an average multiplicative factor on Fq{q, n) is a reasonable approximation 
of "^iFi^q, k), which, if exhibited in detail, would be overwhelmingly complicated without 
rendering significant elucidation or accuracy to justify showing them. 

The results of our calculation of all six terms shown in Fig. 2 will be exhibited below. 
The same procedure can be applied to the determination of hyperon spectra, but will not 
be pursued here. 



IV. TWO-PARTICLE CORRELATION 



A minijet refers to a cluster of particles generated by a semihard parton that emerges 
from the medium. Those particles must be correlated since each of those hadrons in the 
cluster must consist of at least one shower parton generated by the same semihard par- 
ton. The correlation among those multiple shower partons results in correlation among the 
hadrons, even though the hadronization process may involve thermal partons. High-p^ jets 
are routinely studied by jet algorithms, which are, however, ineffective for minijets. Since we 
do not pursue the issue of angular correlation in this paper, that having been done already 



in Ref. 



13| . we focus here on the correlation in the pt variables. The rj and variables 



of the hadrons under consideration are in close proximity, since the correlations are either 
within one jet or in overlapping adjacent jets. Only averages over rj and (p are calculated. 



Let us define the Pearson's covariance, as used in 
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14 



17|, 



= [pi(l)p!(2)]V2 ' 2) = /'^(l, 2) - Pi(l)pi(2) , (9) 

where 

Pi 1 = — 1— , P2l,2 = — — . (10) 

Pidpi PidpiP2dp2 

We shall calculate P2(Pit)P2t) """^ ^^^d pp correlations, but discuss mainly the former. 
We shall use pt and pa (instead of pi and P2) to denote the momenta of the two hadrons {t 
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for trigger, a for associated particle, although the two particles are treated on equal footing 
in p2(l, 2)) in order to avoid notational confusion with the parton momenta pi already used 
in Sees. II and III and in Fig. 1. 

The background subtraction in Eq. is automatically taken into account if we consider 
only the non-factorizable terms in p2{l,2). There are six such non-factorizable terms, which 
are shown schematically in Fig. 3. They are denoted as: (a) (TS)(TS), (b) (TS)(SS)^-^, 
(c) (TS)(SS)2^ (d) (SS)i^(SS)i^ (e) [SSy^iSSf^, and (f) (SS)2J(SS)2J'. The corresponding 
equations for the correlated distributions are given in Appendix A. 




(TS)(TS) (TS)(SS)i^' (TS)(SS)2J 

(a) (b) (c) 




FIG. 3: (Color online) Diagrams for two-pion correlation with same line-types as in Fig. 1. 

Consider first (TS)(TS) correlation. The two shower partons must come from the same 
semihard parton in order for them to be correlated. S'^'^{pi,p2, defined in Eq. I-(A12) 
describes the distribution of those two shower partons pi and p2, integrated over all q > pi + 
P2- Note that their correlation is due primarily to momentum constraints, although charge 
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correlation (or more generally quark-type correlation) cannot be excluded. For simplicity, 
we ignore the constraints arising from quark types, since their recombination with other 
uncorrelated partons (specially thermal ones) partially neutralizes the effect. The parton 
correlation is transmitted to the two pions by the coalescence of the shower partons with 
the thermal partons. 

In (TS)(SS)^-' the correlation arises from the shower parton in (TS) being emitted by 
the same semihard parton that fragments to the other pion. It is obvious that since all 
the diagrams in Fig. 3 are connected, the two pions produced are not factorizable, and 
thus correlated. The equations in Appendix A constitute the six terms of 6*2(1,2). The 
demominator of P2{1,2) in Eq. must include all the factorizable terms also, involving 
not only (TT)j, but also (TS)i(TS)2, etc., where the two shower partons are generated by 
two independent semihard partons. 

Since it has been shown in js] that the proton spectrum in the region px < 5 GeV/c is 
dominated by TTS recombination, and that the p/n ratio decreases with increasing pt above 
the peak at around ~ 3 GeV/c, to study the two-particle correlation between protons it 
is sufficient to investigate only the term (TTS) (TTS). The diagram for that is shown in Fig. 
4 and the equation for it is given in Appendix A. 



P 



P 




(TTS) (TTS) 

FIG. 4: (Color online) Diagram for two-proton correlation with TTS component only for each 
proton. 
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V. PARAMETERS 



There are two parameters adjusted to fit the ALICE data on vr, K,p, A spectra for pt < 
5 GeV/c 7]. They are the inverse slope T of the thermal parton distributions and the 
parameter k in Eq. Their values are js] 

T = 0.38 GeV, k = 2.6. (11) 

Using these parameters we can calculate the spectra at higher including the multi-minijet 
contributions. Although the data for the identified hadrons are not available for pt > 5 
GeV/c, the charged hadron spectrum has been measured up to 19 GeV/c. Assuming that 
S+ can be represented by A, the sum of what we can calculate, it + K + p + A, can be a 
good representation of the charged distribution, as we have demonstrated for p-r < 5 GeV/c 
in js]. Thus we proceed and carry out the calculation for the full pt range. We discover, 
however, that the resultant charged spectrum is lower than the data at high pt- There is a 
reason for that, as we now describe. 



In extending to _^igh px we have fixed k to be constant at the value given in Eq. (fTT]) that 
is determined in |8i]. It has the physical meaning that is the fraction of the semihard 
parton's momentum k that is retained by the parton as it reaches the medium surface, 
since k = nq. Keeping k, fixed implies that the effect of the medium on the energy loss 
is independent of the parton momentum. That property is, however, not consistent with 
the data on the nuclear modification factor Raa at LHC. For 0-5% centrality in Pb-Pb 
collisions at 2.76 TeV ALICE has shown that Raa increases from 0.14 at = 6 GeV/c to 



0.35 at Pt = 19 GeV/c [18|, as is expected from most energy-loss models. It means that 
momentum degradation as a fraction of the initial momentum decreases with increasing 
parton momentum. Thus it is necessary for us to consider a g-dependent k, for which we 
use the form 

where kq and ni are constrained by k = 2.6 at low g ^ 10 GeV/c. We find by fitting the 
charged hadron spectrum that their suitable values are 

kq = 3, Ki = 0.0018 {GeV/c)-\ (13) 
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In the next section we show the results of our calculation based on these parameters and 
compare them to the data. 

Although the overall pt dependence is obtained by varying ko,i, the relative magnitudes of 
the various components are not independently adjustable. Thus we learn about the nature of 
the hadronization processes that we cannot otherwise. Moreover, the correlation properties 
can then be determined without further arbitrariness. 

We note that since K{q) cannot be less than 1, there is an upper limit of for which Eq. 
( !T2|) can be used. It corresponds to g = 33 GeV/c. Any contribution from q larger than 
that value would be very small so the invalidity of Eq. ([T^ for very large q has been ignored 
in our calculation. 




(GeV/c) 

FIG. 5: (Color online) Charged particle distribution in calculated as the sum of four components 
for Pb-Pb collisions at 2.76 TeV. The data are from Ref. [gl. 



VI. RESULTS ON HADRON SPECTRA 



We first show the charged hadron distribution which we identi 
and S. Each of the identified hadron spectra is calculated as in 



as the sum of n,K,p 
but with all multijet 
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contributions included here and with pt extended to 16 GeV/c. The k, parameter used is 
given in Eqs. f|T2|) and ffT3|) . The result is shown in Fig. 5, in which the black solid line is the 
sum that fits the data [gI very well. The four hadronic components are shown separately by 
different lines. Note that for px > 5 GeV/c the pion spectrum dominates over the others, 
but for ~ 3 GeV/c all four have nearly the same magnitude, resulting in the total to be 
noticeably larger than each individually in that region. 
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FIG. 6: (Color online) Meson and baryon spectra at 2.76 TeV with various thermal and shower 
components shown for (a) n, (b) K, (c) p, and (d) S. The data are from Ref. Q] where A distribution 
is shown in (d). 

For mesons and baryons we show them separately in the four panels of Fig. 6. The data 



are from Ref. I?] for low pt, there being none at pt > 5 GeV/c. In Fig. 6(d) the data are for 
A which we regard as indicative of charged hyperon S. The black solid lines in each of the 
four panels are our predictions of the hadronic spectra for 5 < < 15 GeV/c. One of the 
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main points in Fig. 6 is the display of 2-jet contributions to the various spectra. Note that 
for TT in (a) and K in (b) the (SS)^-^ component is always less than (SS)^-', but for p in (c) 
and S in (d) the (SSS)^-' components are almost as large as (SSS)^-' in the 2 < pj' < 6 GeV/c 
region. The ratios (SSS)^-'/(SSS)^-' in the p and S spectra are shown in Fig. 7. At p-p ^ 2 
GeV/c the peaks reach as high as ~ 1. Although T(SS)^'' is not greater than T(SS)"'^-', they 
are approximately equal in the low pt region. (SSS)^-' is small enough to be neglected at all 
Pt- The conclusion is that 2-jet contributions to baryon production, though not large, can 
make quantitative difference in comparing theoretical results to the data. 




6 8 10 12 14 16 

(GeV/c) 



FIG. 7: The ratio of SSS recombination in two jets to that in one jet for p (solid) and S (dashed) 
production at 2.76 TeV. 

To find a phenomenological description of minijets, we recall the experimental way dis- 
cussed in Refs. 
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-l2l[|. where transverse rapidity is introduced, defined by 
yt = lia[{mT + PT)/m^] , 



(14) 



where m.,^ is used in the denominator and in m^, by definition, whether or not the hadron 
refers to a pion or a proton. In those references a peak in yt is found that is attributed 
to minijets. It can be modeled by a Gaussian distribution centered at yt ^ 2.7. Since 
dyt/dpT = l/riiT, the distribuitons dNh/prdpr that we have determined can be readily 
transformed to dNh/ytdyt- Leaving out the TT and TTT components, we show them in 
Fig. 8 for TT and p. We see that various components involving S indeed show peaks in yt 
and that the total of all those components add up to what appear as Gaussians centered 
at yt ~ 2.2 for pion and ~ 2.8 for proton. The rise and fall of the distributions in yt 
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FIG. 8: (Color online) Distributions in transverse rapidity yt for tt and p production at 2.76 TeV. 

are partly due to the definition of yt- The Jacobian relating dNh/ytdyt to dNh/prdpr is 
J = niTPT/yt, which vanishes as pt 0, so even a pure exponential in dN/pTdpr would 
show a bump in dN/ytdyt. The peaks we see in Fig. 8 are, however, due more significantly 
to the suppression of low-pj shower partons. As discussed in [s], the energy loss of semihard 
partons while traversing the medium leads to the enhancement of thermal partons, so the 
partition between thermal and shower partons at low pi is model dependent. The cut-off 
factor 72(^1) given in Eq. I-(A5) marks the end of the shower partons at low pi, resulting in 
the suppression of all TS, SS, TTS, TSS and SSS components at low pr, and hence low yt. 

It is evident from Fig. 8(a) that (SS)^-' is small compared to TS in tt production, and that 
(SS)^-' is even smaller. For p proton production we see from Fig. 8(b) that TTS is dominant, 
but TSS^-' is of comparable magnitude compared to TSS^-' , though both are small compared 
to TTS. The components SSS"-' are significantly smaller in the yt region shown, though 
dominant for > 6 GeV/c. 

The yt distributions in Fig. 8 have the virtue of exhibiting clearly the relative magnitudes 
of the different components in the peak region, more so than in px- The relationship between 
Pt and yt is shown in Fig. 9 where the shaded regions correspond to where the peaks of tt and 
p in yt occur. It is then clear that those peaks are at pt ~ 0.6 and 1.1 GeV/c, respectively, 
although they cannot be easily identified in the plots of dN^/prdpr- 

It is informative to compare the magnitudes of the minijet contributions to tt and p with 
the contributions from thermal partons only. The former are summarized by sum-S depicted 
by solid blue lines in Figs. 8(a) and (b) and are reproduced by the same line type in Fig. 10, 
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FIG. 9: (Color online) pT vs yt showing pion peak in wider (blue) region and proton peak in 
narrower (red) region. 




FIG. 10: (Color online) yt distributions at 2.76 TeV for (a) vr and (b) p, showing sum of all S 
components (blue solid) compared to pure thermal components (purple dashed). Thicker (black) 
solid lines are the totals. 

in which the latter are represented by the purple dashed lines [TT for pion without resonance 
in (a) and TTT for proton in (b)]. The total of all components are the thick black lines in 
Fig. 10. Note that TT distribution in yt has a dip at yt = even though its pt distribution 
is exponential. Except for yt < 1.3, sum-S is larger than TT and TTT, although for the 
latter the proton RF suppresses all components. The implication of this result is that the 
effect of minijets on hadronic spectra is dominant over essentially all pt- 



In the two-component model 19|, |20] a separation of the hadronic spectra into soft and 
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hard components is carried out by demanding that the soft component has a specific form, 
viz., a Levy distribution in mx- What remains after the subtraction of that is the hard 
component, which is characterized (apart from a centrahty factor) by a universal function 
HNNiVt) that is Gaussian in exhibiting a peak at yt ~ 2.7. The region around the 
peak is regarded as being dominated by minijets. The separation is done at the hadronic 
level, working with observed data on pion and proton production at RHIC. What we have 
done in this article is to work at the partonic level, finding many ways in which the shower 
partons can contribute to hadrons, viz., in terms of TS, SS, etc. When the resultant hadronic 
distributions without the soft TT components are plotted in y^, we find the peak originally 
found in Refs. 19|,|20|]. Thus there is phenomenological agreement on what may be identified 



as minijet, although the avenues of approaches to that common ground are quite different. In 
particular, we assume the thermal parton distribution to be exponential that can give a good 
description of both pions and protons at low pt when combined with sum-S. We calculate 
the shower parton distribution but use a cut-off to keep the unreliable part at very low pi 
from exceeding the thermal distribution. Our hadronization scheme treats fragmentation as 
a part of recombination (SS and SSS), and includes cross terms TS, TTS and TSS explicitly. 
Furthermore, by working with the transverse momenta Pi of the partons we have additivity 
that yields the hadronic pt directly. It is a simple property that relates partons to hadrons, 
but is lost in yt- It should finally be noted that, since the TT and TTT components in 
Fig. 10 show peaks in yt without any contributions from shower partons, any claim on the 
existence of minijets cannot be made convincingly without exhibiting correlation, a subject 
to be considered in Sec. VIII. 



VII. HADRONIC SPECTRA AT 5.5 TEV 



At y/sNN = 5.5 TeV the hard parton distribution fi{k) at creation is changed from what 



it is at lower energy and is parametrized in Ref. [22|]. Other parameters, notably kq and ki, 
may differ from those given in Eq. (fT3!) . but without data they cannot be determined. To 
see the relative importance of the different components, it is useful to calculate them for 5.5 
TeV without assuming any specific changes of those parameters. The results on the vr and 
p spectra are shown in Figs. 11(a) and (b). In comparing Fig. 11(a) to Fig. 6(a) for pion 
production we note that while the behavior at low pt is unchanged due to the assumption 
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FIG. 11: (Color online) px distributions for (a) vr and (b) p at 5.5 TeV. 




FIG. 12: The ratios at 5.5 TeV: (a) (SS)2V(SS)iJ' for pion production and (b) (SSS)"V(SSS)iJ for 
proton production with n = 2 (solid) and n = 3 (dashed). 

of unaltered thermal partons, the distribution at the high-p^ end increases by roughly an 
order of magnitude. Furthermore, (SS)^-' become almost as large as TS for Pt ^ 6 GeV/c. 

A more visibly apparent change is seen in proton production by comparing Fig. 11(b) to 
Fig. 6(c). The (SSS)^-' component increases by almost two orders of magnitude at pt ~ 16 
GeV/c. Even at ~ 2 GeV/c its value at 5.5 TeV is enhanced over that at 2.76 TeV by 
more than a factor of 10. One sees in Fig. 11(b) that multijet recombination is important at 
all Pt- To compare the various multi-minijet contributions, we show in Fig. 12(a) the ratio 
(SS)2V(SS)i^' for pion production and in Fig. 12(b) (SSS)2V(SSS)i^' and (SSS)3V(SSS)i^' for 
proton production. Note that the proton peak for (SSS)^-'/(SSS)^-' at 5.5 TeV in Fig. 12(b) 
is more than twice higher compared to that at 2.76 TeV in Fig. 7. 
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FIG. 13: p/tt ratio at 2.76 and 5.5 TeV. 
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FIG. 14: (Color online) Shower parton contributions to (a) pion and (b) proton distributions in yt 
at 5.5 TeV. 

Since both pion and proton distributions at large pt are higher at 5.5 TeV compared to 
those at 2.76 TeV, the p/vr ratio is not significantly changed; both are shown in Fig. 13. The 
peak is shifted slightly lower to around px — '^ GeV/c. 

The distributions in yt at 5.5 TeV have the general structure as at 2.76 TeV, as shown 
in Figs. 14(a) and (b) for pion and proton, respectively. The peaks occur at around the 
same values of yt, however, the magnitudes are much higher — nearly triple. Notice that 
(SS)2J' is now as high as (SS)i^ while (TSS)^^' exceeds (TSS)i^'. At both 2.76 and 5.5 TeV, 
the sums of all shower parton contributions (shown by the solid blue lines in Figs. 8 and 14) 
are significantly higher than the (SS)"-' and (SSS)"-^ contributions that do not involve the 
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FIG. 15: (Color online) Sums of S contributions compared to (a) TT and (b) TTT distributions 
in yt at 5.5 TeV. 

thermal partons. The differences between sum-S and pure T components are shown in Fig. 
15(a) for TT and (b) for p; evidently, at 5.5 TeV sum-S is now more than 4 times larger than 
TT for vr at the peak, and as much as 6 times TTT for p. Clearly, the minijet contributions 
dominate over pure thermal on the one hand and over pure fragmentation on the other. 



VIII. RESULTS ON TWO-PARTICLE CORRELATIONS 

For Tx-Ti correlation we use Eqs. (A1)-(A6) to calculate dN^^T, / ptPadptdpa- There is no 
TT term in either the trigger or the associated particle because it is factorizable. All six 
terms included are shown in Fig. 3 and are clearly non-factorizable. The result on the 
correlation function P2^{ptiPa) defined in Eq. (I9l) is shown in Fig. 16(a) for collisions at 2.76 
TeV, exhibiting clearly that the correlation is important only when pt and pa are in the 1 
GeV/c region. The most important contribution to that correlation in the peak region comes 
from the (TS)(TS) term in Fig. 3(a) that involves two shower partons produced by a single 
semihard parton. That cluster of partons with low transverse momenta is what has been 
referred to as minijet, and the correlation is among the components of that cluster, but for 
correlation between pions there is a gain in momentum due to recombination of the correlated 
shower partons with thermal partons in the vicinity. The properties near the correlation 
peak can best be shown in a plot of P2^{yti,yt2) in Fig- 16(b). The peak is located at 

= yt2 ~ 2.2. This result is consistent with the findings at RHIC in Au-Au collisions at 200 
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FIG. 16: (Color online) tt-tt correlation at 2.76 TeV on (a) pt-Pa and (b) yt^-yt2 




FIG. 17: (Color online) (a) Ratios of components in tt-tt correlation: (TS)(SS)^''/(TS)(TS) and 
(SS)"'^-'(SS)"'^-'7(TS)(TS) on pt-pa, (b) P2{yti,yt2) ^oi a pion from (SS)^-' correlated with another pion 
in terms shown in Fig. 3 (c), (e) and (f). 



GeV for unlike-sign correlation on the same side Uj . To show the suppressed contributions 
from the other components compared to (TS)(TS), Fig. 17(a) exhibits (SS)^-'(SS)^-' and 
(TS)(SS)^-' relative to (TS)(TS). They are obviously insignificant for the values of pt and Pa 
shown. For (SS)^-' contribution to P2{yti,yt2) shown in Fig. 17(b), they are also very small 
compared to that shown in Fig. 16(b), which is dominated by (TS)(TS). 

The correlation between two protons has even more components. We show in Fig. 18 only 
the dominant component (TTS)(TTS), as depicted in Fig. 4; here the same peak is shown 
in plots on (a) pt-Pa and (b) yt^-yt^. Compared to tttt correlation in Fig. 16, the magnitude 
of pp correlation is an order of magnitude smaller. The peak in yt^-yt2 is shifted to slightly 
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FIG. 18: (Color online) (TTS)(TTS) contribution at 2.76 TeV to p-p correlation on pt-Pa and (b) 

yti-yt2- 

higher value at ~ 2.8. 

The above results are our best theoretical evidence for minijets. There are, of course, 
jets produced at all higher px. However, since inclusive two-particle distribution involves 
integration over the parton momenta qi, higher-p^ jets would not show up unless cuts in pt 
are made to select those jets. When that is done, one should be able to obtain peaks in the 
autocorrelation on the angular-difference variables \rii — 772! and |0i — 02|, observed in the 



data 



11 



3, 



23| and calculated in 



A summary of our findings in this section is that the peaking of two-particle correlations 
in TTTT and pp distributions provides indisputable evidence for minijets whose correlated 
shower partons are responsible for the phenomenon, as observed in [1JL|. Without particle 
identification, the correlation is due mainly to the process depicted in Fig. 3(a). Two-jet 
recombination is of negligible effect in correlations. 



IX. CONCLUSION 



We have calculated the hadronic spectra for meson and baryon production in Pb-Pb 
collisions at LHC for p^ in the range 0-16 GeV/c. Of particular concern in this work is the 
investigation of the extent to which multi-jet recombination is important in that pt range, 
since jet density is high at LHC. We have found that at 2.76 TeV the (SS)^-' component 
in meson production makes negligible contribution at all pt compared to other components 
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which are primarily TT and TS at low pt, and then (SS)^-' for pt > Q GeV/c. For baryon 
production (SSS)^-' is comparable to (SSS)^-' at ~ 2 GeV/c, so jet fragmentation alone 
in the conventional sense is not reliable at such low pt- Recombination involving T is more 
important for pt < Q GeV/c; in particular, (TSS)"-' components with n = 1 and 2 have 
comparable magnitudes, either one being larger than (SSS)"-'. Since that is the low-p^ 
region where the B/M ratio is of the order of 1, the two-jet contributions to the inclusive 
spectra should not be ignored. The effect of 2-j recombination to the two-hadron correlation 
is, however, negligible. 

It is not surprising that at 5.5 TeV the effect of multi-minijet contribution becomes 
more important. For proton production (TSS)^-' is as large as (TSS)^-' for p^ < 5 GeV/c 
and (SSS)^-' is comparable or larger than (SSS)^-' for nearly all px- At ~ 6 GeV/c all 
components except TTT and (SSS)^-' are approximately similar in strength. Thus we predict 
that any calculation without 2-jet recombination would not be able to reproduce the data. 

We have shown by studying the single-particle distributions in yt that there are peaks at 
i/t ~ 2.2 and 2.8 for tt and p, respectively, verifying what has been observed experimentally by 



STAR 



19142 1|. Because of the definition in yt at low yt, a Gaussian distribution in yt does not 



unambiguously indicate the existence of minijets. However, when two-particle correlation 
exhibits a broad peak in yt-yt distribution, then the inference of minijets is inevitable. What 
we have shown is that the peak in i/t-yt receives its dominant contribution from (TS)(TS) in 
the case of tttt correlation, and from (TTS)(TTS) in pp correlation. That is our explanation 
of the data on that peak observed by STAR The corresponding value of pt where the 
peak is located is around 1 GeV/c. Thus the minijets that give rise to the shower partons 
have a significant effect on low-p^ physics through thermal- shower recombination. 

The dominance of the role played by minijets in the low-pr region puts a new light on the 
subject of soft physics, since the conventional treatment by hydrodynamics does not take 
minijets into account explicitly. It raises the question on whether the non-flow component 
can ever be ignored at high collision energies, especially at LHC. So far we have examined 
only 0-averaged distribution in pt for central collisions. The problem that lies ahead is 
clearly the study of azimuthal anisotropy for non-central collisions. It is there that we must 
confront the LHC data on f„ j25|, and clarify the roles of minijets versus the fluctuations 
of initial-state configurations, both of which appear to have similar effects on the final-state 
hadronic observables. 
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Appendix A: Two-particle Correlations 

We summarize in this Appendix the formulas for the various terms contributing to tt-tt 
correlation, followed by one term for p-p correlation. With background subtraction defined 
for (72(1,2) in Eq. ([9]), it is sufficient to list only the non-factorizable terms involving S in 
various combinations. The functions k), T(pi), Si{x), S'^{pi,k) and S'^'^{pi,p2, n) used 
in following equations have been defined in Eqs. I-(38), I-(5), l-(Bl), I-(A4) and I-(A12), 
respectively. 
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